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Abstract 

We consider a class of toy models where a spatially flat universe is filled 
with a perfect fluid. The dynamics is found exactly for all these models. In one 
family, the perfect fluid is of the phantom type and we find that the universe is 
first contracting and then expanding while the dynamics is always accelerated. 
In a second family, the universe is first in an accelerated expansion stage, then 
in a decelerated expansion stage until it reaches a turning point after which it 
contracts in a decelerated way (increasing contraction rate) followed by another 
accelerated stage (decreasing contraction rate). We also consider the possibility 
to embed this perfect fluid in a realistic cosmology. The first family cannot be 
viable in a conventional big bang universe and requires a rebound in the very 
early universe. The second family is viable in the range < 1+to^o < 0.09 for 
a spatially closed universe with a curvature satisfying current bounds. Though 
many of the models in this family cannot be distinguished today from a universe 
dominated by a cosmological constant, the present accelerated expansion is 
transient and these universes will reach a turning point in the future before 
entering a contraction phase. 

PACS Numbers: 98.80.-k, 98.80.Es, 95.36,+x 
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1 Introduction 



Since the discovery of the universe present accelerated expansion rate pQ, our view 
of expanding universes has been durably affected. There has been a lot of interest 
in models producing the accelerated expansion [2] as we see it at the present time. 
The simplest models introduce dark energy (DE) in the form of a comoving perfect 
fluid with a sufficiently negative pressure today. A very simple and appealing exam- 
ple is provided by a cosmological constant A which agrees well with existing data. 
Remarkably, this corresponds to a constant energy density and a constant pressure 
satisfying p A + = or = — 1- Leaving aside the conceptual problems raised 
by the existence of such a tiny cosmological constant, it remains to be seen whether 
more accurate observations will force us into a paradigm shift. In any case it remains 
an important issue to look for other viable alternatives where the energy density and 
pressure will have a different physical interpretation and behaviour. In particular it 
may well be that no DE component is needed at all if the present accelerated expansion 
is caused by modifications of gravity on cosmic scales. A miscroscopic interpretation 
of the DE energy density and pressure is given in models like quintessence [3], a 
time-dependent homogeneous scalar field minimally coupled to gravity. Through and 
Einsteinian interpretation of the modified Friedmann equations, a microscopic de- 
scription of an effective DE energy density and pressure can be given as well for some 
models going beyond General Relativity, (see e.g.[l] in the context of f(R) models) 
where gravity is modified on cosmic scales. This radical departure from conventional 
expansion histories has also brought interest for DE of the phantom type (wde < — 1 ) 
[5] and for the future of our universe [6]. As we will see below, our work touch upon 
both issues. 

Here we introduce a class of phenomenological perfect fluids and we study the 
corresponding universe dynamics. The crucial property is that the energy density 
of this perfect fluid can vanish in the course of expansion or contraction, depending 
on the specific model. We consider toy models where spatially flat universes contain 
only this perfect fluid and we solve for their dynamics. As expected this gives rise to 
very peculiar dynamics. In one family the perfect fluid is forever of the phantom type 
and it will tend to a cosmological constant in the asymptotic future after passing 
through a bounce. In the second family the universe expands till a turning point 
after which it starts contracting. We also consider the possibility to include such 
a component in a realistic cosmology. We will find that viable universes can exist 
with a radical departure from standard cosmological scenarios either in the past, for 
the perfect fluid of the phantom type, or in the future for the second family. The 
departure in the universe past in the first case comes from a conceptual problem and 
requires probably to go beyond General Relativity in the early universe. However for 
the second family a universe with a peculiar future can be obtained in the framework 
of General Relativity. 
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2 The model 



We consider a Friedmann-Lemaitre- Robertson- Walker (FLRW) spatially flat universe 
filled with a comoving pefect fluid of energy density p and pressure p. We have the 
following equations 

& = *-fp (i) 

H = -4ttG (p + p), (2) 

The time evolution of the perfect fluid 

p= -3H(p + p) = -3H(l + w)p , (3) 

is contained in (jTJ),([2J) and we have introduced the equation of state (EoS) parameter 
w = p/p. Let us start with a perfect fluid satisfying 

p = constant . (4) 

It is easy to derive from ([I]),® 

d(p a 3 ) = -p da 3 = -d(p a 3 ) , (5) 

which leads to 

.3 



or 



Setting p = —A, we can write 



p a 3 = -pa 3 - B , (6) 
B 

P = -p-->0. 7 

a 4 

. B 

p=A--. 8 

a 6 

It is easy to understand why both terms can appear in (jSJ): from Q, the first term 
leads to p and p being constant while the second term corresponds to p — 0. The 
important point here is that B is an arbitrary integration constant. It can vanish and 
while it does not, it can have arbitrary sign, the sign is not specified by (HI). 

The case B > corresponds to a phantom equation of state, w < —1. This 
will always be the case, there is no crossing of the phantom divide line. We have in 
particular from (j5J) 

H > . (9) 

It is not possible in this case to interpret (jSJ) as a cosmological plus some dust because 
the additional dustlike component would have a negative energy density due to the 
"wrong" sign of B. It is the sum of both terms in (jSj) that gives a physically sound 
energy density satisfying p > 0. 

Solving for the scale factor a as a function of t we obtain 



as cosh 3 



(10) 
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Here a B is the minimal value of the scale factor reached at time t B where the Hubble 
parameter H vanishes. Equation (JSJ) can be recast as 



P 



A 



1 - 



a 



Though the universe we consider is spatially flat, we get a bouncing universe with 
contraction at times t < t B and expansion at times t > t B . We have defined 



:i2i 



which represents the value of H in the asymptotic future. The Hubble parameter H 
in function of t is readily obtained from ([TO]) , viz. 



H = Hnr, tanh 



- H^it — t B ) 



(13) 



from which we see again that H tends to for t — > oo. To summarize, our universe 
is always in a phantom stage, it is first contracting then expanding and tends to a 
pure de Sitter universe with H = for t — > oo. 

It is interesting to compare ffTDl) and f[T3T) with the solutions obtained in a matter 
dominated universe with a cosmological constant A 



n 



m,0 



a sinh3 ( - H^t 



H 



coth ( ^ H^t 



(14) 
(15) 



^. In this case we have of course a "Big Bang" singularity with H — > oo 



with Hi 

as t — > 0. It is the past which is radically different for both universes. 
We can look for more general behaviours 



P 



A - — 

a n 



AB > . 



(16) 



We have in particular A — p > 0. It is easy to obtain the corresponding pressure p 



n 



P 



(17) 



We consider first the case n > 0. We check that a phantom type equation of state 
is obtained for n > 0, a fact already clear by inspection of ([IB"]) . Hence a universe 
filled with a perfect fluid that evolves according to (JTBjl undergoes only accelerated 
expansion. We can generalize the expressions (TT0|) . ( II 3p . viz. 



as cosh™ 



i^oo tanh 



n 



— H^t — t B ) 

n . . 

— Hoc\t ~ tB) 



(18) 
(19) 
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Actually these expressions apply to negative n as well. It is interesting to consider 
also the case n < as well. Of course, now we have no phantom type equation of 
state and H < The universe is a de Sitter space in the infinite past and we can call 
its (constant) Hubble parameter H^. It will reach some turning point in the future 
with expansion followed by a contraction stage. We have now (for n < 0) 



H 



ob cosh 
if™ tanh 



n 



n 



Hoo(t — t B ) 



Hoo(t — t B ) 



(20) 
(21) 



It is seen that a_ n = a n 1 from which we obtain, in accordance with eqs. ({TB"]) and 

(EUD 

(22) 



H- n = —H n 



For n < 0, H < so that the expansion is not always accelerated. Actually the 
evolution is now very different. The universe is first expanding and then it starts 
contracting at time t B . In the asymptotic past, t — > — oo, the universe tends to 
an expanding de Sitter space while it tends to a contracting de Sitter space in the 
asymptotic future, t — > oo. During some (symmetric) interval around t B one has a < 
similarly to a spatially closed universe going over from expansion to contraction. At 
some time t a < t B , the evolution goes over from accelerated expansion to decelerated 
expansion on the interval t a = t B — At < t < t B with At = t B — t a . This is followed 
by decelerated contraction on the interval t B < t < t B + At. For t > t B + At the 
universe undergoes accelerated contraction. The interval At is easily found, viz. 



sufiV 



I 71 1 

T 



(23) 



It is seen that the interval At depends on |n|. Hence for n < 0, the universe evolution 
can be suumarized as follows: 



a > 
a < 

with At found from ( 123]) . and further 

a > 
d < 



\t - t B \ > At 
\t - t B \ < At 



t < t B 
t>t B . 



Let us turn our attention to the equation of state (EoS) parameter w 
straightforward to derive 



n (A 



w 



- 1 



(24) 
(25) 



(26) 
(27) 

It is 

(28) 



3 VP 

which shows that w diverges whenever H = 0. This has to be the case as we show 
now. 
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Assuming a contracting universe (n > 0) starts from some regular initial conditions 
at time t\ and reaches a bounce at time the following equality must hold from 

mm 



' 5 V 



dt (1 + w) = -oo . (29) 



ti 



Clearly, this equality cannot be satisfied for tg < oo if w is regular. Assuming an 
expanding universe with n < starts from some regular initial conditions at time t\ 
and reaches a turning point at time tg, we get analogously 



t B 

dt(l + w)=oo, (30) 



leading to a similar conclusion. 



3 Embedding in a realistic cosmology 

A further interesting issue concerns the embedding of our perfect fluid in a realistic 
cosmology. This is what will be investigated now. It is natural to assign to our perfect 
fluid the role of dark energy responsible for the late-time accelerated expansion of the 
universe. Hence our fluid will be called below dark energy with energy density p^E 
and pressure Pde- For our discussion it is convenient to introduce the notation 



B = B , 



(31) 



so that f lT6|) can be recast in the form 



p = A-B - 

where Bq has dimensions of energy with Bo = A — po- It is easy to show 



(32) 



a 



Bo 
A 



(33) 



When n > 0, as corresponds to a bounce with — > for Bq — > 0. When n < 0, as 
corresponds to a turning point with a# — > oo for Bo — > 0. The Friedmann equations 
yield now 



H 2 

a 
a 



8vrG 
47rG 



Pm + Pr + A - B [ — 

a 



[Pm + 3p r - 2PDE ~ Tl(A - p DE )) 



(34) 



(35) 



where we have allowed for spatially non-flat universes. We finally note that n = 
corresponds to a cosmological constant, only the difference A — Bo matters with 
8ttG (A — Bo) = A and this case needs no further consideration. 
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3.1 n>0 



When n > 0, the universe will tend in the future to a de Sitter space. The dark energy- 
sector is of the phantom type, it tends asymptotically to a cosmological constant 
from below (wee < —1) which is not possible for usual quintessence (a minimally 
coupled scalar field) . We have curiously a "growing" rather than "decaying" effective 
cosmological constant. 

The past of the universe must be changed compared to the standard big bang 
scenario. While our perfect fluid has no influence of the universe dynamics in the 
past - actually it becommes dynamically significant only at z < 1 - the universe 
is not allowed to cross in the past the value a = a B if we insist on the condition 
Pde > 0. On the other, if we want a universe where the main features of the hot big 
bang universe is retained, B has to be small enough so that 



which requires a fine tuning of the model in addition to the fine tuning of A (the 
cosmic coincidence). Viability of such a universe requires further that the universe 
has a bounce at a > a B when p DE is still positive. This requires to go beyond 
General Relativity in the very early universe (e.g. [7]. In any case, due to (1561) . such 
a universe is undistinguishable observationally from a universe with a cosmological 
constant. We believe the main merit of this model is that it provides an interesting 
conceptual problem. 



When n < 0, the primordial universe can be governed by equations fl34l) . fl35|) with a 
standard big bang cosmology. The universe future is changed with a turning point 
after which the universe will be constracting. 

The danger with such a universe is that a B could lie in our past which is forbidden 
as further expansion would force pde to be negative and the universe could not reach 
the present day expansion stage. The condition B < A is obviously sufficient to 
guarantee that a B lies in the future, a B > a . Taking ^ < 1 will push a B (far) 
in the future. In such a model dark energy is represented by a slowly "decaying" 
cosmological constant. Actually, even if a B 3> a , once the point a = a B is reached, 
further expansion - which is unavoidable because of the presence of dustlike matter 
(and radiation) however tiny their energy densities - would inevitably imply pe>e < 0. 
This is reminiscent of the situation encountered with a Big-Brake singularity in the 
future [8] if one is willing to include dustlike matter. 

The only way to avoid this is to have a closed universe. Let us call a x the value of 
the scale factor at which the dustlike term and the curvature term cancel each other 
in the right hand side of the Friedmann equation. We then have 




(36) 



3.2 n < 




(37) 
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Inclusion of radiation is possible but it will have negligible influence. To ensure that 
the turning point is reached before pde becomes negative, we must impose a\ < as- 
This can be expressed as a lower bound on jf^o | in function of the model perameters 

|^fc,o| > ^m,o (^"^j " ' fity 

In a realistic cosmology, —0.0065 < Qk,o < 0.0012 (95% C.L.) with central value 
^fc,o = —0.0027 [9j. So there is preference for a slightly spatially closed universe 
which is exactly what we need. Condition fl38|) can obviously be satisfied for a large 
range of parameter values satisfying 

< (0.026) |n| . (39) 
where we have defined Qbo = - ^ iL ~ an d &ao = and we have taken f2 m0 — 0.25. 

' Pcr,0 ' Pcr,0 ' 

Violation of fl38|) for some set of parameter values implies that the observational upper 
bound on spatial curvature is not large enough to ensure pde > in the future. 
Inequality ( 1391) can be recast into 

i^2_ < \ (40) 

^de,o ~ (0.026)-l"l - 1 v ' 

For example for \n\ = 2 we have Qb,o < 0.00067 £Ide,o, f° r \ n \ — 0.5, one obtains 
Qb,o ^ 0.192 Qde,o- For |n| < 0.5, Qb,o can be substantial while it becomes negiligible 
for |n| > 2. 

On the other hand, the expansion must be accelerated at the present epoch which 
yields 

2Vl£)Efi ~ ^m,0 ^ 1-25 

iZ s,o < rn ~ t~t > ^ 4i J 

|n| |n| 

Condition (T4"ll corresponds to wde,o < ~1^deo (neglecting here f2fc )0 ). We know 
however from the observations that wde,o is significantly smaller. Hence a much 
stronger condition comes from the smallness of 1 + wee,o- Indeed we have 

"D£,0 \ n \ 

If we take 1 + wee,o ^ 0.1 and Qde,o — 0.75, we obtain Qb,o ^ nfjF- As expected, 
this is much stronger than ( l4"Tl) . 

Combining condition (H01 with (1421) . we get a new inequality 

77* 1 

~L + w DE0 < — — n (43) 

,u ~ 3 (0.026)-l"l - 1 V ' 

It is interesting that (}4"3~|) implies 

1 + wm,o < 0.09 . (44) 
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Hence this model is able to cover the range of values relevant for 1 + wde,o > 0. As 
the observational bound on |f4,o| will be more stringent, the range of allowed values 
for 1 + Wde,o > will decrease. 

Let us consider the regime |n| > 2. In that case we have Qb,o < 6 x 10~ 4 £Ia,o, 
Qb,o < 6 x 10~ 4 Qde,o and Qa,o — ^de,o- We also get in this case from (H3"j) that 1 + 
wde,o < 4x 10~ 4 . These models cannot be distinguished from a cosmological constant 
model with Qa,o — ^a,o- Actually for models with |n| > 0.95, the quantity 1 + wde,o 
is strongly constrained with 1 + wde,q ^ 0.01. For large |n| values both quantities 
and j^ 2 ^ are much smaller than one. It is possible to have substantially higher 
values for 1 + wde,o when |n| —> 0. For |n| ~ 0.5, the quantity 1 + %g i0 is allowed to 
be inside the range < 1 + wde,o ^ 0.03. For lower nonzero values |n| ~ 0, 1 + wde,o 
can be as large as 0.09. In this range of small \n\ values, „ B '° is allowed to be much 

larger than 1 with f2 B0 ~ Qa,o (and f2 B0 < Qa,o, ^a,o — ^b,o = ^de,o ~ 0.75). The 
quantity fi^o can be large because |n| is very small, and therefore the time evolution 
of Pde is slow enough so that one can still obtain a# ^> ao and a# > 

If observations show evidence for wde,o significantly higher than —1 but still 
satisfying wde,o ^ 0.09, then only the low \n\ values remain viable. In all these 
universes the late-time accelerated expansion is a transient phenomenon. It is followed 
by decelerated expansion and eventual contraction. 



4 Summary 

We have considered spatially flat toy models universes filled with a particular perfect 
fluid and we have solved exactly the corresponding dynamics. In one family of models 
(n > 0) the perfect fluid is of the phantom type and an ever accelerating dynamics is 
obtained: the universe is first contracting, then passing through a bounce after which 
it is expanding. We have exponential contraction in the infinite past and the universe 
tends to a de Sitter expanding universe in the asymptotic future. In the second 
family of models (n < 0) the universe is exponentially expanding in the asymptotic 
past and exponentially contracting in the asymptotic future. At some time (t a ) the 
expansion becomes decelerated, the scale factor then reaches a maximum followed 
by a contracting phase, first decelerating then accelerating (decreasing contraction 
rate). We have also considered the possibily to embed our perfect fluid in a realistic 
cosmology in which case it can play the role of dark energy. When (n > 0) we have 
the strange situation that while the corresponding energy density is negligible already 
at low redshifts, a standard big bang singularity would force the energy density of our 
perfect fluid to become negative. We conjecture that this could be avoided by having 
a rebound before this happens, which could be the case if we go beyond General 
Relativity in the very early universe. On the other hand, at the present time this 
universe cannot be distinguished from a A dominated universe and it tends to a de 
Sitter (expanding) universe in the future with 1 + wde < 0. 

As for the other family of solutions (n < 0), a slightly spatially closed universe 
is required in order to avoid pde to become negative in our future. Interestingly, 
using the current observational bound on jO^ol we find that all universes with < 
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1 + wde,o J$ 0.09 can be accomodated by this model. The maximal values 1 + 
WDEfl — 0.09 are obtained for small values ^ \n\ <C 1. For |n| > 2, we have 
essentially a A dominated universe with Qb,o *C ^a,o — ^a,o and 1 + wde,o <C 
1. The strong inequality fi B <C Qa,o can be seen as an additional fine tuning. 
These models certainly do not pretend to get rid of any of the conceptual problems 
of the A model, what is interesting is that they can provide a radically different 
future. Because the effective cosmological constant is slowly decaying, the accelerated 
expansion is a transient phenomenon and the present expansion stage will be followed 
by a contracting phase even though this universe cannot be distinguished today from 
a A dominated universe. The future of our universe may depend crucially on the 
precise form of the dark energy sector even if it cannot be distinguished today from 
a cosmological constant, and our model provides just such a example. 
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